An Electromagnetic Signature of Supermassive Black Hole Binaries 
That Enter Their Gravitational- Wave Induced Inspiral 
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Mergers of gas-rich galaxies lead to black hole binaries that coalesce as a result of dynamical 
friction on the ambient gas. Once the binary tightens to < 10 3 Schwarzschild radii, its merger is 
driven by the emission of gravitational waves (GWs). We show that this transition occurs generically 
at orbital periods of ~1-10 years and an orbital velocity v of a few thousand km s _1 , with a very 
weak dependence on the supply rate of gas (v oc M 1//s ). Therefore, as binaries enter their GW- 
dominated inspiral, they inevitably induce large periodic shifts in the broad emission lines of any 
associated quasar(s). The probability of finding a binary in tighter configurations scales as v~ 8 owing 
to their much shorter lifetimes. Narrow-band monitoring of the broad emission lines of quasars on 
timescales of months to decades can set a lower limit on the expected rate of GW sources for LISA. 

PACS numbers: 04.25.dg,98.62. Js,98.54.Aj,98.62.Ra 



Introduction. Recent advances in general relativistic 
simulations of the final coalescence phase of black hole 
(BH) binaries through gravitational wave (GW) emis- 
sion sparked new interest in related electromagnetic 
(EM) signals 0. New observational searches were initi- 
ated for wide BH binaries in their pre-merger phase or 
single quasars which are displaced from their host galaxy 
spatially or spectroscopically due to GW recoil after the 
merger [3j. In addition to testing the recent numerical 
predictions, such searches can be used to calibrate the 
expected rate of GW sources for the proposed LISA mis- 
sion as well as to reduce the error box on the source 
localization after LISA is launched. The next generation 
of simulations will attempt to incorporate gas dynamics 
and radiative transfer in order to calculate the EM coun- 
terparts to GW sources If these EM signatures are 
sufficiently unusual relative to typical quasar variability, 
one could search for them well before LISA is launched. 
Since BH binaries spend more time at larger separations, 
wider binaries should be more abundant in quasar sur- 
veys. In this paper we study a generic EM signature of 
the widest quasar binaries whose dynamics is dominated 
by GW emission. 

The formation of BH binaries is a natural by prod- 
uct of galaxy mergers (which are generic to the hierar- 
chical build-up of galaxies in the standard cosmological 
model) , since almost all galaxies in the local Universe are 
observed to harbor a BH at their center [6|]. The tidal 
torques generated during a merger extract angular mo- 
mentum from any associated cold gas and concentrate 
the gas near the center of the merger remnant, where its 
accretion onto a BH results in a bright quasar [t| . 

Binary Parameters and Statistics. We consider two 
black holes with masses Mi and Mi in a circular orbit of 
radius a about each other 0] . Their respective distances 
from the center of mass are a% = (fj,/Mi)a (i — 1,2), 
where fx = M\MijM and M = Mi + M 2 . We define the 
parameter £ = 4/i/M, which equals unity if Mi = M2 



and is smaller otherwise. The orbital period is given by, 



P = 2ir(GM/a 3 



-V = 1.72 yr a% 2 M~^, 



(1) 



where, a 16 = (a/10 16 cm) and M s = (M/10 8 M Q ). The 
angular momentum of the binary can be expressed in 
terms of the absolute values of the velocities of its mem- 
bers vi and t)2 as J = £j=i,2-Mji>jaj = fMva, where the 
relative orbital speed is 



vi + v 2 = (27TO/P) 



1.15 x 10 4 km s _1 M 8 1/2 o 16 1/2 . 

(2) 

We focus our discussion on mergers of gas-rich galax- 
ies since those are most likely to result in bright quasars 
which are easy to observe. Hydrodynamic simulations 
of such mergers indicate that the associated BH binaries 
inspiral steadily as a result of dynamical friction on the 
ambient gas [9(. In contrast, BH binaries in gas-poor 
mergers, which shrink by scattering stars, may stall once 
they reach a separation a < lpc if the supply time of new 
stars which cross their orbit exceeds the Hubble time (re- 
sulting in the so-called "final parsec problem" [13]). In 
gas-rich mergers, the rate of inspiral slows down as soon 
as the gas mass interior to the binary orbit is smaller 
than fi and the enclosed gas mass is no longer sufficient 
for carrying away the entire orbital angular momentum 
of the binary, J [ll| . Subsequently, momentum conserva- 
tion requires that fresh gas will steadily flow towards the 
binary orbit in order for it to shrink. The binary tight- 
ens by expelling gas out of a region twice as large as its 
orbit (similarly to a "blender" opening a hollow gap [121 ]) 
and by torquing the surrounding disk through spiral arms 
[HElI El- Fresh gas re-enters the region of the binary 
as a result of turbulent transport of angular momentum 
in the surrounding disk. Since the expelled gas carries 
a specific angular momentum of ~ va, the coalescence 
time of the binary is inversely proportional to the sup- 
ply rate of fresh gas into the binary region. In a steady 
state, the mass supply rate of gas that extracts angular 
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momentum from the binary, M, is proportional to the 
accretion rate of the surrounding gas disk. Given that a 
fraction of the mass that enters the central gap accretes 
onto the BHs and fuels quasar activity H, 17 1, it is ap- 
propriate to express M in Eddington units M. = M /Me, 
corresponding to the accretion rate required to power the 
limiting Eddington luminosity with a radiative efficiency 
of 10%, M E = 2.3M yr^Mg. We then find, 



i gas w ( J/Mva) 



Li/M = 1.1 x 10 7 yr (M~ 



(3) 



For a steady M., the binary spends equal amounts of time 
per log a until GWs start to dominate its loss of angular 
momentum. 

The coalescence timescale due to GW emission is given 
by 0, 



5 c 5 a 4 



2.53 x 10 5 yr 



(4) 



256 G 3 M 2 ^ j~ ^ M 3- 

By setting tow = tgas we can solve for the orbital speed, 
period, and separation at which GWs take over, 

v GW = 7.2 x 10 3 km s" 1 C 1/4 (MM a ) 1/8 ; (5) 

F GW = 7.1 yr C 3/4 M 8 5/8 ^- 3 / 8 ; (6) 

a GW = 8.3 x 1CT 3 pc C 1/2 Ml /4 M- 1/4 . (7) 

For a binary redshift z, the observed period is (1+z)Pqw- 
Remarkably, the orbital speed at which GWs take over is 
very weakly dependent on the supply rate of gas, i>gw x 
M 1 / 8 . It generically corresponds to an orbital separation 
of order ~ 10 3 Schwarzschild radii (2GM/c 2 ). 

In an equal-mass binary, each BH moves relative to 
the center-of-mass at a speed of |i>gw- Coincidentally, 
this orbital speed happens to be similar to the veloc- 
ity width of the broad emission lines (BEL) of quasars 
[19j |. This fortunate coincidence implies that the time- 
dependent gravitational potential and UV illumination 
of the binary members will introduce variability into the 
spectroscopic shape of any associated BEL. It is believed 
that BEL are emitted by numerous compact clouds which 
cover about a tenth of the sky of the quasar d<Jj. Even 
if the BEL clouds originate in the outer edge of the ac- 
cretion disk around one of the BHs, their dynamics will 
be affected by the combined gravitational force and UV 
illumination from both BHs. 

We therefore robustly predict that when BH binaries 
from mergers of comparable-mass, gas-rich galaxies start 
to inspiral by GW emission, they would show variability 
with a period of months to decades (for M(l + z) 8 / 5 ~ 
10 6 -10 9 Mq) in any associated BEL, across a range of 
Doppler velocities comparable to the BEL spectral width. 

The notion that BH binaries would display periodic 
shifts in their emission lines was recognized a long time 
ago [2()1 |. So far, intensive monitoring programs for a 
limited number of active galactic nuclei revealed spectro- 
scopic variability but no reproducible evidence for a pe- 
riodic signal [2l|. The new insight added by Eqs. (07]) 



is that BH binaries should generically start their GW- 
induced inspiral at orbital periods that are observation- 
ally accessible, and at orbital velocities that are compa- 
rable to the BEL width. This would naturally lead to a 
double-peaked line profile which can be observed to vary 
periodically. In contrast, the line profile is expected to 
remain singly peaked [38| at the lower orbital velocities 
associated with wider binary separations (where the su- 
perposed BEL of the two BHs are only slightly offset in 
velocity relative to each other), or for binaries that are 
much more compact than their common BEL region (and 
hence act gravitationally as a single BH). 

The periodic Doppler shift of the BEL is expected to be 
smaller for unequal mass binaries. If each BH controls the 
spectroscopic centroid of its own BEL feature, the BEL 
associated with M, will be shifted by v; = (/J,/Mi)v. The 
relative abundance of BH binaries with different mass 
ratios can be modeled theoretically [22I ]. 

The probability of observing periodic BEL shifts due 
to a binary transitioning to a GW-dominated inspiral, 
?gw, depends on the unknown fraction of the associ- 
ated quasar lifetime which overlaps with this final phase 
of binary coalescence. In case accretion is suppressed 
by the opening of an inner cavity in the disk [l2j , it is 
possible that no BEL emission is observable early in the 
GW-dominated inspiral phase. However, recent numeri- 
cal simulations 13 , Hi , 15 , 3, 13 indicate that accretion 



should continue during this phase through the formation 
of two secondary accretion disks around the individual 
BHs. If quasar emission persists throughout the coa- 
lescence process at a constant M. due to the abundant 
supply of cold galactic gas, then one is likely to find as 
many binaries with separation oqw ~ 10~ 2 pc as with a 
separation of tens of pc (for which the enclosed gas mass 
is comparable to li). In that case, the prospects for dis- 
covering a quasar binary with a periodic BEL shift are 
good, since an example for a tight binary at a projected 
separation of ~ 7.3 pc has already been identified Q 
and three sub-pc candidates have been reported recently 
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The probability of finding binaries deeper in the GW- 
dominated regime, V cx <gw, diminishes rapidly at in- 
creasing orbital speeds, with V = Pgw^/^gw) -8 - 

Discussion. State-of-the-art numerical simulations 
15L 17 1 focused so far on wide binary separations where 



gravitational torques dominate the interaction of the bi- 
nary with the surrounding gas, but did not attempt to 
describe the late, viscosity-driven, evolutionary phase of 
interest here. This late inspiral phase is mediated by 
angular momentum transport in the circumbinary disk 
which delivers fresh gas from large distances into the bi- 
nary region. Since quasar accretion disks on scales of 
~ 10 4 Schwarzschild radii are radiation-pressure dom- 
inated, a proper numerical treatment of the turbulent 
transport of angular momentum requires a radiation- 
magneto-hydrodynamics code. Suitable numerical sim- 
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ulations were only applied to a single BH so far 25j , but 
could be extended with a global code [26| to the binary 
accretion problem. 

Given the limitations of existing simulations, one is 
restricted to modeling the inflow of gas into the binary 
orbit through analytic approximations. A recent analytic 
treatment of viscous transport within a circumbinary ac- 
cretion disk [lO ] ignored the existence of individual com- 
pact accretion disks around the two BHs and adopted 
a simplified prescription for a-viscosity and binary mi- 
gration. The formation of two compact accretion disks 
around the BHs (that is necessary for producing BEL) 
is expected to occur through gas inflow from the inner 
edge of the circumbinary disk 27j , and the moving inner 
disks could affect the gas flow around the binary both 
hydrodynamically (l6j and radiatively (especially if they 
radiate near the Eddington limit [lltl- Bearing these 
simplifications in mind, Fig. 3 in Ref. [ll| is in excellent 
agreement with Eq. ([5|), confirming that the GW-driven 
phase is approached at a characteristic orbital velocity of 
~ 10 4 km s _1 with only a weak dependence on the ac- 
cretion model details. The generic results of Eqs. ([5][7]) 
naturally lead to double-peaked line profiles which vary 
periodically on timescales of months to decades. 

Future monitoring of the broad emission lines of 
quasars can set a lower limit on the expected rate of GW 
sources for LISA. A multi-year search dedicated to peri- 
odic variability of BEL in a large sample of quasars has 
not been performed as of yet. Preliminary studies have 
already made use of spectroscopic data from SDSS [28| 
and SDSS-III [29J. Any related discovery would motivate 
the implementation of narrow-band filters in photomet- 
ric variability surveys, culminating with the future LSST 
[30I ] . A precursor for this approach may be realized in 
the near future. 

The Palomar Transient Factory [3l[ (PTF) plans to ac- 
quire Ha filters at 6564A which are 80-100Awide. The 
periodic Doppler shift of the BEL in quasar binaries that 
enter their GW-dominated inspiral corresponds to a frac- 
tional wavelength change of (vgw/c) ~ 1% or ~ 66A. It 
would be possible to notice such a shift as an observed 
line slides in and out of the above mentioned filters. The 
sky surface density of quasars brighter than an i-band 
magnitude of 19, for which the Mg II, C IV or C III 
lines are redshifted to 6564 A, is - 0.2 per square de- 
gree on the sky. This results in 3 quasars per pointing 
if two filters are used. The planned survey will moni- 
tor up to several hundred fields containing up to ~ 10 3 
quasars with a cadence of 1 month, and measure relative 
i-band fluxes to a precision of 0.3% (3 milli-magnitudes). 
It is also possible to design a dedicated survey that will 
cover 40% of the sky (a few thousand quasars) with a 
1-year cadence. Separately from this program, a survey 
is underway to study variability (reverberation mapping) 
of quasars brighter than an r-band magnitude of 21 us- 
ing six 200A wide filters of the COMBO-17 survey ^ 



on the Chandra Deep Field South, using the ESO 2.2 m 
telescope at La Silla [331 ] . 

Spectroscopic surveys are much more labor intensive. 
Previous spectroscopic searches for BEL variability did 
not report any detection of periodicity for several tens of 
quasars that were monitored over many years 34| or in a 
sample of a few thousand quasars that were sampled only 
a few times over less than a year 3|| 3f|, but identified 
examples of longer-term variability without a measured 
period [2p| . 

Since quasar BEL are known to vary for other reasons 
211 ] , it is crucial to identify periodicity before classifying 
a candidate source as real. In confirmed cases, measure- 
ments of the period and velocity shift can be used to 
solve for M and a up to an ambiguity regarding the in- 
clination angle of the orbital plane relative to the line of 
sight. This ambiguity may be removed through reverber- 
ation mapping [19] , |36j , since a measurement of the time 
delays between temporal variations in the UV continuum 
originating near the BHs and the two response times of 
their BEL, would provide an independent estimate of Mi 
and Mi (and hence £). Even before periodic variability 
is observed, reverberation studies can be used to distin- 
guish the double line profile of a binary (for which the 
two lines would respond differently to variations in the 
UV continuum of the two quasars 37J) from that of a 



single disk emitter (for which the two lines will vary si- 
multaneously). In addition, the quasar luminosity could 
be used to estimate M based on related hydrodynamic 
simulations [H El EE EE Ezl ■ With all this informa- 
tion, observers will be in a position to examine whether 
i gas is indeed of order £gw m these BH binaries. 
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